Permanents in linear optical networks 
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We develop an abstract look at linear optical networks from the viewpoint of combinatorics and 
permanents. In particular we show that calculation of matrix elements of unitarily transformed 
photonic multi-mode states is intimately linked to the computation of permanents. An implication 
of this remarkable fact is that all calculations that are based on evaluating matrix elements are 
generically computationally hard. Moreover, quantum mechanics provides simpler derivations of 
certain matrix analysis results which we exemplify by showing that the permanent of any unitary 
matrix takes its values across the unit disk in the complex plane. 
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I. INTRODUCTION 

Quantum computing promises to be able to perform 
certain computational tasks such as factoring of large 
numbers [ll exponentially faster than by using classical 
computing. Such statements can be made when look- 
ing into computational complexity classes of these al- 
gorithms. Recently, a wealth of proposals for physical 
implementation of such algorithms has emerged. The 
particular one we have in mind is based on qubit en- 
coding in photons which are subjected to passive lin- 
ear transformation in beam splitter networks. These 
so-called quantum-gate engineering protocols make use 
of measurement-induced nonlinearities in which effective 
nonlinear evolutions are achieved by conditional partial 
measurements in an extended Hilbert space augmented 
by auxiliary modes [J . The probabilistic nature of these 
schemes makes it necessary to optimize linear-optical net- 
works with respect to their probability of success. Hence, 
apart from the complexity of the algorithm that is sup- 
posed to run on the network there is another complexity 
associated with the optimization problem. 

It seems to have been first noted in |3| (and some- 
what more discussed in 4]) that optimization of those 
networks with respect to their probability of success is 
closely linked to maximizing a permanent under certain 
constraints. Permanents have found their way into quan- 
tum physics by the work of Caianiello [5|, Ig who showed 
that generally all expectation values of bosonic field op- 
erators can be written as a permanent (or a hafnian, a 
matrix function related to the permanent which was ac- 
tually introduced by him) of a certain coefficient matrix. 
At that time it seemed to be little more than a notation 
although he was eventually able to sum up Dyson series 
in field theory without the use of Feynman graphs. 

In this article we will show that indeed the natural way 
of writing unitary transformations of multi-mode Fock 
states is by means of permanents of matrices closely as- 



sociated with the unitary matrix of the state transfor- 
mation. This immediately implies that the calculation 
of matrix elements of unitarily transformed multi-mode 
states (and all subsequent calculations based on it) is 
computationally expensive and cannot be performed in 
polynomial time (polynomial in the number of modes). 

The article is organized as follows. Section ^1 is used 
to introduce some definitions and notations. The main 
result of this paper, the unitary transformation of multi- 
mode states, is then presented in Sec. lIIII As an example 
for the application of these formulas, in Sec. IIVI we de- 
scribe how the entanglement power of linear-optical net- 
works can be calculated. The paper finishes with some 
concluding remarks in Sec. 



II. SOME DEFINITIONS AND NOTATIONS 

For the sake of definiteness, we introduce some nota- 
tion which will be used throughout this article. To begin 
with, let us specify what one commonly understands by 
the permanent of a matrix. We will concentrate only on 
square matrices since this is what we will be using later. 

Def. : The permanent of an [n x n)-matrix A is a gen- 
eralized matrix function defined by 



per A = ^ W^ia, 



(1) 



a£S„ i=l 
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where Sn is the group of permutations. For a fixed per- 

n 

mutation a Cz Sn, the product Y[ ^iui is called a diago- 

i=i 
nal. The permanent is thus the sum over all diagonals of 

the matrix A Q- 

It is easy to see that the permanent is similar to the 
determinant of a matrix, but lacks the signs related to the 
order of the permutation. In fact, these signs cause the 
permanent to be a very special matrix function in that 
all usual computational rules known for the determinant 
fail. For example, no rules known from linear algebra can 
be used to simplify the computation of the permanent. 
Especially, rules such as det OSO = det O det S det O 



= det S for orthogonal matrices O do not hold. The 
only known simplification for permanents can be achieved 
for permutation matrices P and diagonal matrices D in 
which case we have 



pcrPAD = perPper AperD 



(2) 



In general, the computation of the permanent has to be 
done without possible simplification in which case there 
are n! terms in the sum (^. In fact, there exists a result 
by Marcus and Mine @ that asserts that there exists no 
linear transformation on n x n-matrices with n > 3 such 
that the permanent of these matrices can be converted 
into a determinant. 

It will prove to be convenient to introduce some more 
notation. Let A[fci, . . . , km\h, ■ ■ ■ , Im] be the {m x ni)- 
matrix whose matrix elements are those of the original 
matrix A with row indices fci , . . . , km and column indices 
Zi, . . . , Im- For example. 



Mki,k2,k3\h,h,l3] = 



Afc,i, 



Afc2i2 




(3) 



The object A[(l'"i , 2™^ ...) 1(1"! , 2"^ ...)] denotes a ma- 
trix whose entries are taken from the matrix A and whose 
row index i occurs exactly rm times and whose column 
index j occurs exactly rij times, for example. 



A[{1\2\S')\{V\2\3')] 



Ai2 Ai2 Ai3 

A22 A22 A23 I . (4) 
A32 A32 A33 



Furthermore, let Gn,N be the set of all non-decreasing 
integer sequences a;. 



G„ 



,N 



:= {uj : 1 < uji < ■ ■ ■ < ujn < N} . (5) 



That is, each sequence u) £ Gn.w has length n and its 
entries take integer values up to N. 



III. UNITARY TRANSFORMATION OF 
MULTI-MODE FOCK STATES 

With the above definitions it is now rather straightfor- 
ward to write down the expression how a bosonic (pho- 
tonic) Fock state \ni,n2, ■ ■ ■ ,71^) transforms under the 
action of a unitary network specified by an operator U 
or a unitary matrix A, respectively. Traditionally, this 
transformation is written by using the transformation 
rule for photonic amplitude operators |^] 



A^at 



(6) 



U\ni,n2,.-. ,njv) 



uMi 



N 

E 



Afc 




By using the multinomial expansion theorem, this can be 
further rewritten as 

U\ni,n2,---,nN) ^ ^ ^ 



N 



n "»j! 

N N 

n (A,aa]j"«i . . . n (A,„A.a]„)««^« |0)®^ . (8) 
Note that the n^ form a n x n-matrix whose row sums 

N 

are the photon numbers rii, hence ^ riij — rii. More- 

over, the row sums (as well as the column sums to ensure 
unitarity, hence conservation of the total photon num- 
ber) add up to the total photon number n which makes 
the matrix {uij/n} (doubly) stochastic. Denoting the 

N 



column sums of {f^^j} as ^ 
second line in Eq. © as 



riii = m 



■J' 



we can write the 



/ N 

n 



1/2 



N 



mi\ 



Y[ Al\^'\\mi,m2,...,mN). (9) 



,k,l=l 



N 



Note that Y[ ^7k' ^^ ^ product with exactly n factors 

where each row index i occurs rrii times and where each 
column index j occurs Uj times. Therefore, it must be a 
diagonal (according to the definition in Sec.|nJ of the ma- 
trix A[(l™i,...,iV™«)|(l"i,...,iV"")]. If we take the 
permanent of this matrix, we see that out of all possi- 
ble permutations of column indices, Y[j ''^r of those per- 
mutations are identical. Analogously, there are Oi"^*' 
ways of distributing the row indices. Hence, not all di- 
agonals are different from each other. However, only 
(rii^^iOdli "•jO/dli 7 "yO terms actually lead to the 
same diagonal. With the knowledge of these multiplici- 
ties it follows from Eq. (jHJ that 



{mi, 1712,.. .,mN\U\ni,n2,. .. ,nN) 




-1/2 



-1/2 



n^ 



perA[n'|n] (10) 



where 



n = (i"\2"%...,iv""), 

O' = (1™S2'"%...,7V™") 



(11) 
(12) 



Equation (|10|l is proof for the intimate relation between 
unitary transformation of multi-mode Fock states and 
permanents of matrices associated with the unitary ma- 
trix A. 

An immediate consequence of Eq. \iQ\ is that 



per A 



®N , 



l\U\l 



,®N 



(13) 



with the effect that the permanent of a unitary matrix 
takes its value across the unit disk in the complex plane, 
I per A I < 1, since the expectation value to find the state 
|1)®^ after performing a unitary transform of a simi- 
lar state is actually a probability amplitude. Although 
this result is usually derived from the Marcus-Newman 
theorem pJJ |, the probability interpretation of quantum 
mechanics provides a much simpler derivation. 

The description H10|l of matrix elements can be 
straightforwardly generalized to the full unitary trans- 
formation as 



C/|ni,n2,. .. ,nAr) = ]^?ij! 



-1/2 



(14) 



1 



V J_^ perA[a;|r?]|mi(a;),m2(a;),...,mjv(a;)) 

in which the mi{ijj) are the multiplicities of the occur- 
rence of the value i in the non-decreasing integer sequence 
u) and ^.{u>) = Y\^ mi{u))\. 

Analogously, we obtain partial matrix elements (or 

rather projections) with X]i=2 '^j < n as 



112,... 


,mN\U\ni,n2,. . 


.,nN) = 


(1 


N 


\ ^^'^ i "" \ 


-1/2 


N 


[W 


11-.! 


peTA[n"\n] 


n-> ]m, 


i=l 


/ \,=i ) 




i=2 



with 



o" = (i"-i:,=2'"j-^2'^ 



,iV™") 



(16) 



Equations H1U|) and (|14|) constitute the main result of 
this paper. Either equation implies that quantum-state 
transformations can be written in terms of permanents 
of matrices associated with the unitary matrix A. Al- 
though nothing has changed in the state transformation 
itself, this result nevertheless shows that the computa- 
tional complexity of computing matrix elements in the 
Fock basis and functions depending on it is the same as 
computing permanents. The intrinsic interest in this re- 
sult stems from the fact that it provides an argument for 
the affiliation of any type of optimization algorithm for 
passive linear optical networks or any maximization or 
averaging procedure over unitary operation in the Fock 
basis to the same complexity class, namely that of com- 
puting a permanent. The following section shall provide 
an example for an averaging procedure over unitary op- 
erator that can entirely be written as a function of per- 
manents. 



IV. ENTANGLEMENT POWER OF UNITARY 
NETWORKS 

Over the recent years, entanglement and its quantifi- 
cation has played a major role in the discussions on the 



foundations of quantum mechanics. The basic principle 
lies in the distinction between separable and inseparable 
(entangled) states. Commonly, a bipartite state is called 
separable if its density matrix can be decomposed into 
a convex sum of tensor product states of the respective 
subsystems. On the contrary, an entangled state cannot 
be written in this form. Apart from separability criteria 
that can tell to which of the above classes a given quan- 
tum state belongs, there exist entropic or distance-based 
measures to quantify the amount of entanglement. We 
will not give any details here since they can be found in 
the vast literature on this subject (for a recent review, see 
e.g. [ij), we merely note that for pure bipartite states 
all of these entanglement measures collapse to the von 
Neumann entropy of the reduced density matrix. 

A somewhat related question is to ask how much en- 
tanglement a quantum operation can generate when ap- 
plied to some initial bipartite state |12|. Among all pos- 
sible questions of this type we will concentrate only on 
the problem how much entanglement will be generated 
on average from a pure product state as opposed to the 
maximally available entanglement [13. We will hence- 
forth call a functional Pe{U) the entanglement power of 
the unitary operator U with respect to the entanglement 
measure E if 

PEiU) ^ J dgi\i;i))dgi\i;2)) E{U\i^,,^2)) , (17) 

where dg{\ipi)) denotes the group measure with respect 
to the state \ipi). In what follows, we will restrict our- 
selves to the linear entropy as the entanglement measure 
of choice, simply because it is easiest to compute. The 
linear entropy, defined by 



L(L/|^i,V2)) = 2(l-Tr^?) 



(18) 



where gi = TT2{U\i^i,fp2){'^i,fp2\U^), although not a 
proper entanglement measure in the strict sense, serves 
as an upper bound to the distillable entanglement. 



A. Qubit initial states 

To begin with, let us consider the simplest example in 
which two qubit states impinge on a beam splitter with 
transmittivity T £ C We will not use the language of 
permanents here but rather perform a straightforward 
calculation because that appears to be simpler (we leave 
the calculation involving permanents to the example in 
the next section). Let us define the input states as 



IV'2) 



co|0) ■ 
do\0) 



ci|l) 
-dill) 



(19) 
(20) 



where |coP + |cip == jdoP + \di\^ = 1, and let the unitary 
matrix associated with the beam splitter be given by 



A = 



T R 

-R* T* 



(21) 



The usual replacement rules for the bosonic amplitude 
operators, a^ i— > A"^a^, immediately imply that 



t/|^i,V2)=codo|00) 

+ icodiR + cidoT)\10) + (codiT* 

+cidi \V2TR\20) - V2T*R*\02) 



(22) 



cidoR*)\01) 

^i\Tf~\Rf)\n) 



The cocfBcicnt in the last term in Eq. H22() is propor- 
tional to the permanent of the beam splitter matrix 
per A = |Tp — |i?p, as one expects from the general 
theory presented above (it is just the probability ampli- 
tude of finding the state |11) from a |11) input state). 
The reduced density matrix of subsystem 1 is then 



91 = 



V2cidiT*R*\0) 



h.c. (23) 

[{codiT* - cidoR*)\0) + cididTp - |i?n|l)] ® h.c. 
codolO) + {codiR + cidoT)\l) + V2cidi\2)] (g) h.c. . 



The remaining steps are to compute the linear entropy 
and to integrate over the group measures with respect to 
the dynamical groups associated with the initial states 
\ipi)- The crucial observation at this point is that the 
group integration removes all phase dependencies from 
the linear entropy, an effect that can be attributed to 
Schur's Lemma. The remaining amplitude dependence 
(either on \T\ or on |_R|) can be uniquely rewritten in 
terms of the permanent, since jTp = (1 -t- per A)/2 and 
|i?P — (1 — perA)/2. A straightforward calculation then 
leads to the result that 



Pl{U) 



64 



[l - (perA)2] [13 + 9(perA)2] (24) 



which only depends on the permanent of the beam split- 
ter matrix. Obviously, if the beam splitter acts only as 
a phase shifter (in which case per A — 1) or as a two- 
sided mirror (per A = — 1) the input beams do not mix 
and no entanglement is created whatsoever. The maxi- 
mum of Eq. (|24|l is reached when per A — 0, i.e. when 
the beam splitter is balanced. The entanglement power 
at this point is 39/64 and monotonically decreasing with 
increasing |perA|. 



B. Higher-dimensional initial states 

Let us now look at the situation in which the initial 
states that impinge on a single beam splitter are A^ + 1- 
dimensional, hence can be written as 



N 



N 



l^l) = Yl Cnil^^l) , 1-02) = Yl ^"2l'^2) , (25) 



ni— n2— 

such that the combined state reads 

N 



|'01,V'2)= Y '^ntdnjni,n2) 

^1,^2—0 



(26) 



Now we refrain from using the straightforward calcula- 
tion since that becomes messy very quickly. Instead, we 
use the permanent language developed in the previous 
section. Applying the transformation rule (|14|l to this 
state yields 



U\iJl,ij2) 



N 



Cn-i 0:n 



ni,n2— 



v: 



71l!7l2! 



v-^ perAfa;|(l"i,2"=)], , , , ,,,^^, 



ajGG„^ + „2,2 



^toi(u;)!to2('^)! 



where we used the notation (1"^ , 2"^ ) to indicate the mul- 
tiplicities of the occurrence of the indices 1 and 2 in the 
columns of the matrix A[a;|(l"i, 2"^)], as described in 
Sec. ^ The important object to compute is now the 
linear entropy. It is straightforward to see that one can 
write the trace over the square of the reduced density 
matrix as 

Tvgl = ^(fci,A;2|f7|Vi,V'2>(^i,V'2|L^t|A:i,fc2) 

ki ,k2 

{i;i,HU^k2){k2mi,^2). (28) 

This shows that we need all ingredients from the previous 
section, that is, both matrix elements of a unitary opera- 
tor and projections or partial matrix elements. It is now 
easy to see that the matrix elements (fci, k2\U\ipi, 1P2) can 
be written as 



N 



{h,k2\U\^Pl,ij2}= J2 



^m ^n2 



^/ ki\k2lni\n2 



perA[(l'=S2'=^)|(l"S2"^)],5„,+„„fe,+fc, (29) 

In complete analogy, we can cast the matrix elements 
{A,i'2\U^\k2) {k2\U\^Pi,^P2) in the form 



(^I>2|f7t|fc2)(fc2|f7|^i,^2) = 



N N 

E E 



mJ^2=On,7^=o ^2!("i + ^2 ~ fe) ! V^i !m2 !ni !n2 ! 

perA*[(l"l+"^-'^^2'=^)|(l"^2"^)] 

perA[(l"l+"^-'=^2^^)|(l"^2"^)]5„,+„„„,+„,.(30) 

Thus, Eq. (|28|) can be written in functional form, by using 
Eqs. (|29|l and l|30|l . in terms of permanents associated 
with the beam splitter matrix A. 

Although the number of sums to be performed to com- 
pute the trace of the square of the reduced density ma- 
trix counts up to ten, there are already three Kronecker 
6 functions that reduce the number of sums to seven. 
Written out explicitly, Eq. (|28|l reads 



Tr^? = E 



AT 



N 



AT 



N 

E 



Cpi dp2 C(ji flgj C,„i "rn2 '^"1 ''" 



perA[(l'=S2'=^)|(FS2P^)]perA*[(l''S2'=^)|(l«S2«=) 
perA*[(l"l+"^-'^^2'=^)|(l"^2"^)]perA[(l"l+"^-'=^2'=^)|(l"^2"^)] 

'^Pl +P2 , A;i +/i;2 '^Ql +^2 1 A;i +/i;2 *^"ll +^TT'2 ,ni +712 

I 



(31) 



At this point we need to have a look at the restrictions 
the group integration imposes. The coefficients Ci, i = 
1, . . . , N, are complex and their squared moduli add up to 
unity. Hence, they can be represented by hyperspherical 
coordinates and some phases as 



on absolute values of certain permanents associated with 
the beam splitter matrix A. 



Multi-mode networks 



cos 0i+i IT sin Qj , Q 



N+l 



0. 



(32) 



With the additional restriction that ipo — (because an 
overall phase can be absorbed into the definition of the 
states) we are left with 2N integration variables per ini- 
tial state. Going back to Eq. (|28|l one notes that the 
integration has to be done over four coefficients of the 
form c*c*CkCi (and similarly for the coefficients of the 
second initial state). Averaging over the phase angles 
then yields 



N 



(2^) 



N 



Y[ d,f^ 



j{(Pi+iPj-iPk-ipi) 



(SikSji + Sii6jk){'^ ~ hSij) 



(33) 



Since this type of integration has to be performed for 
both initial states, there are only few sums remaining in 
Eq. I2H1). 

The integration over the angles Qi is done by noting 
that the coordinates on an A^-sphere S'^ induce a surface 
integration measure 



dt.{s'')=(j[sin^-'e,)(j[de) 



(34) 

(35) 
(36) 



where V is the surface area of the A^-sphere. Eventually, 
we are left with 



'ore, we find that 


3 


(A + l)(A + 3) 
1 


(A + l)(A + 3) 



2 3 



{SjkSji + SuSjk){l + ^S,j) 
(A + l)(A + 3) 



(37) 



Inserting this result into Eq. 1)28(1 , we finally obtain some- 
thing which depends only at most quartically on perma- 
nents of the type |per A[(l'"l+'"2-'=^ 2'=2)|(l'"i , 2"'^)] , i.e. 



A final remark should be made about multi-mode net- 
works. In principle, the calculations can be repeated in 
this case as well and the matrix elements be expressed 
in terms of permanents. However, since it is not en- 
tirely clear how to define entanglement measures (or even 
monotones) in a multipartite setup, we are not attempt- 
ing to use any partial results in these cases. Obviously, 
once a measure for multipartite states has been estab- 
lished, it must be expressible in terms of permanents, as 
we have seen. 



V. CONCLUSIONS 

In this article we have shown that calculation of matrix 
elements of unitary operators in the Fock basis are inti- 
mately connected with the computation of permanents. 
This at first surprising result becomes entirely obvious 
if one realizes that the (re-)distribution of photons in a 
linear-optical network is a purely combinatorial problem. 
As a matter of fact, one can think of these networks as 
being represented by a bipartite graph whose adjacency 
matrix has complex entries. 

With the knowledge about the relation between ma- 
trix elements of unitary operators and permanents, it be- 
comes immediately obvious that the computational com- 
plexity of calculating functions that depend on perma- 
nents is the same as that for permanents. Examples 
for such algorithms that appear in quantum informa- 
tion theory include maximizations of success probabilities 
on linear-optical networks Q or averaging procedures in- 
volving unitary operators as in the entanglement power 
of such networks. 

The advantages of relating unitary transformations to 
permanents are twofold. On one hand, unitary transfor- 
mations of multi-mode photon states can be written in a 
rather compact form. On the other hand, we think it may 
be possible of gaining more insights into the theory of per- 
manents of unitary matrices from a quantum-mechanical 
viewpoint. The example we gave that the permanent of a 
unitary matrix takes its value across the unit disk in the 



complex plane shows that quantum-mechanical reason- 
ing can simplify the derivation of such propositions. We 
believe that this combinatorial view on linear-optical net- 
works in quantum information processing can stimulate 
the emergence of more intricate results. 
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